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Abstract Motivated by the possibility that the fun- 
damental "constants" of nature could vary with time, 
this paper considers the long term evolution of white 
dwarf stars under the combined action of proton de- 
cay and variations in the gravitational constant. White 
dwarfs are thus used as a theoretical laboratory to 
study the effects of possible time variations, especially 
their implications for the future history of the universe. 
More specifically, we consider the gravitational con- 
stant G to vary according to the parametric relation 
G = Go(l + t/t±)~ p , where the time scale t* is the 
same order as the proton lifetime tp. We then study 
the long term fate and evolution of white dwarf stars. 
This treatment begins when proton decay dominates 
the stellar luminosity, and ends when the star becomes 
optically thin to its internal radiation. 
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1 Introduction 

One of the overarching but unresolved questions in 
physics is whether or not the fundamental constants 
are truly constant, or if they could vary with time. 
A related question is whether these constants could, 
in principle, have other values in far-away regions of 
space-time (i.e., in effectively other universes). Current 
experiments place limits on the time variations of the 
fundamental constants, e.g., the gravitational constant 
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G and the fine-structure constant (Uzan 2003 and 
references therein) . These constraints can be expressed 
in terms of corresponding time scales, tq = G/G and 
te = as/atE, which are found in the range 10 10 — 10 12 yr 
(Uzan 2003). Although these time scales are safely 
longer than the current age of the universe (to = 13.7 
Gyr; e.g., Spergel et al. 2007), these values are much 
shorter than some time scales that are experimentally 
accessible; as one example, the proton lifetime has a 
measured lower bound of 10 33 yr (Supcr-K 1999). As 
a result, over the vast expanses of time available to the 
universe in the future, time variations in the physical 
constants could change the projected future history of 
the universe. 

This paper addresses this issue by considering pos- 
sible time variations of the gravitational constant and 
its corresponding effects on the long term evolution of 
white dwarf stars. In this context, we are using white 
dwarfs as a theoretical laboratory to consider the effects 
of time- varying G. This choice is justified because white 
dwarfs are among the simplest and hence most well un- 
derstood stellar objects (starting with Chandrasekhar 
1939) and because they play an important role in the 
future history of the universe (Adams & Laughlin 1997, 
hereafter AL97; see also Cikovic 2003, Dyson 1979, Is- 
lam 1977). For example, almost all stars turn into white 
dwarfs after their nuclear burning phase, and hence a 
sizable fraction of the accessible baryonic content of the 
universe is locked up in white dwarfs after stellar evo- 
lution has run its course (for completeness, note that 
a sizable fraction of the baryons remain in the medium 
between galaxies in dusters, e.g. Nagamine & Loeb 
2004). 

To address this issue, we must define both the time 
scales and the functional form of the time variations 
under consideration. Although an enormous variety of 
time variations are possible, we narrow the scope by 
allowing the gravitational constant G to vary according 
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to the parametric relation 

G(t) = G (l + i/**)" p , (1) 

where Go is the current value, and where t* and p are 
parameters. This functional form is perhaps the sim- 
plest type of allowed time variation — the value of G 
reduces to the current value at the current cosmological 
epoch and decreases as a power-law in the long term fu- 
ture. Note that current experimental limits imply that 
r G = G/G > 10 12 yr (e.g., Table IV of Uzan 2003; see 
also Barrow 1996) so that t* S> to- In the long term fu- 
ture, proton decay is one of the most important energy 
sources for white dwarfs. This problem thus has two 
time scales of interest: the proton lifetime tp and the 
time scale t* for variations in G. In the limit t* 3> tp, 
white dwarf evolution closely follows previous results 
obtained using a fixed gravitational constant (AL97, 
Adams ct al. 1998; Dicus et al. 1982). As a result, 
this paper primarily considers the limit t* <C tp, where 
white dwarf evolution is dominated by changes in grav- 
ity, and the case where t* is comparable to tp. 

Although this work will focus on time variations in 
the gravitational constant G, for completeness we note 
that time variations in other physical quantities arc 
possible. For example, the masses of the fundamental 
particles or, equivalcntly, the ratio of electron to pro- 
ton mass fi = m e /mp, could be time dependent (Cal- 
met & Fritzsch 2006). We note that the quark masses, 
which play a role in determining the proton mass, are 
actually the fundamental constants of the underlying 
theory; however, the ratio \i has been experimentally 
constrained (see Uzan 2003 and references therein). 

In this paper we consider the variation of G by itself 
in order to isolate the effects of its possible time evo- 
lution. However, fundamental theories (string theory, 
M theory) suggest that the forces are unified and hence 
may vary in strength in a coupled manner (Uzan 2003 
provides a brief review of the possibilities). Unfortu- 
nately, we do not have a definitive working theory of 
how such coupled time variations are expected to oc- 
cur, and we leave such a more comprehensive treatment 
for future work. 

This paper is organized as follows. In §2, we out- 
line a basic model for white dwarf structure, including 
our treatment of time variations in G and proton de- 
cay. One interesting complication is that proton decay 
leads to a downward cascade for the atomic weights of 
the constituent nuclei; this issue is addressed in some 
detail. In §3 we outline the basic results of the paper, 
including the H-R diagrams for the long term evolution 
of these degenerate stars and their corresponding chem- 
ical evolution. We then conclude in §4 with a summary 
and discussion of our results. 



2 White Dwarf Model 

2.1 Conditions and Assumptions 

Given the extremely long time scales associated with 
baryon decay and possible time variations of the grav- 
itational constant, we make the following assumptions 
about the properties of the white dwarf stars under con- 
sideration here. The end of conventional star formation 
is expected to occur around ~ 10 14 yr after the big bang 
(AL97), and baryon decay will not occur until a cosmic 
age of ~ 10 33 yr (Super-K 1999) or later. As a re- 
sult, we can safely assume that the white dwarfs will 
cool to low temperatures. Over time scale 10 14 — 10 22 
yr they will be kept warm through dark matter cap- 
ture and annihilation (AL97), but afterwards have only 
proton decay as an energy source. Also, any angular 
momentum the white dwarfs may have acquired during 
formation will be lost, leaving a body free of rotational 
deformations. As a result, electron degeneracy pressure 
alone must act to counterbalance gravity and maintain 
the star's structure against gravitational collapse. 

For white dwarf stars in the low temperature limit, 
as considered herein, the physical structure of the star 
decouples from thermal considerations. We can thus 
assume that the equation of state is that appropriate for 
non-relativistic degenerate matter, so that the pressure 
is related to the density by the relation 

P = Kp 5/3 , (2) 

where the constant K is given by 

K=^(3nr /3 (^e)- 5/ \ (3) 

where m e is the electron mass and m is the atomic 
mass unit (Shapiro & Tcukolsky 1983). Here, we take 
/i e = (A/Z), the average ratio of atomic mass to atomic 
number for the entire collection of particles comprising 
the white dwarf star. This ratio will be a function of 
time due to the effects of baryon decay (see £12.31 ) 

The star must be in hydrostatic equilibrium, which is 
described by the Lane-Emden equation, a dimensionless 
differential equation of the form 

1 d 2 d6 _ 

j£ ~ ~ e ' (4) 

where 

r = a£ , 

P = Pc0 n , 



Long Term Evolution of White Dwarf Stars 



3 



0.8 



0.6 



<X> 



0.4 



0.2 



1 2 3 4 5 6 7 

Fig. 1 Dimensionless form of the density profile given by 
the Lane-Emden equation for polytropic index n = 3/2 
(non-relativistic, solid line) and n = 3 (relativistic, dashed 
line). Note that p tx 9 n . 
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Iii this model, the white dwarf star is a polytropic, self 
gravitating gas, with polytropic index n = 3/2. Figure 
Q] shows the density profile of such an idealized white 
dwarf for the non-relativistic case (n = 3/2). For com- 
parison, the relativistic case (n = 3) is also shown. Note 
that white dwarfs are characterized by n = 3 (relativis- 
tic) polytropcs only for masses approaching the Chan- 
drasekhar mass; the vast majority of white dwarfs pro- 
duced by the universe will have lower masses and can 
be described as n = 3/2 polytropes (e.g., Shapiro & 
Teukolsky 1983). 

The Lane-Emden equation provides a satisfactory 
description for the pressure gradient in the interior of 
the system where the pressure is extremely high. How- 
ever, once the star has radiated away all of its residual 
thermal energy, its surface temperature will be nearly 
absolute zero. For the outer layers of a frozen white 
dwarf, the typical inter-particle separation is about one 
Bohr radius, and the constituent particles naturally fall 
into a tightly packed lattice structure. This densely 
packed state for the matter in a frozen white dwarf ef- 
fectively places a positive, non-zero lower bound on the 
density. Here, we take the lower bound on the density 
of this lattice to be a constant value po = 1 g cm -3 . 
At this density the coulomb, self-gravitational, and de- 
generacy energies have comparable magnitudes and the 
electrons are no longer degenerate. 

Early in a white dwarf's lifetime, this non-degenerate 
surface layer is extremely thin, with little of the stellar 



mass residing therein, and it is a reasonable approxi- 
mation to assume that the entire system is degenerate. 
For instance, early in the life of a 1.0 Mq white dwarf, 
the non-degenerate 'crust' is only about 1 km deep. On 
the time scales of this paper, as the star expands due 
to mass loss from proton decay and varying G, it is es- 
sential to follow the boundary between degenerate and 
non-degenerate layers within the white dwarf. 

2.1.1 Mass 

The total remaining mass of the white dwarf at any 
time is given by the equation 



M(t) = M exp[-rt] 



(5) 



where T is the baryon decay rate. For the sake of def- 
initeness, we take T = 10~ 37 yr _1 for this work; the 
issue of proton lifetimes is discussed in $2.'2\ 

Invoking conservation of mass and the definitions 
of equation (J4j) , the mass contribution of the degen- 
erate core enclosed within the degenerate surface, £,deg, 
is given by 



M deg (t) = r) m [p c (t) G(t) - 3 fi e (t) 



-51 V2 



(0) 



where 
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'deg\ 



An alternative way of calculating the total mass is to 
substitute £i and |#i(£i)| for ^deg and |6^ eff |, respec- 
tively, for the expression r\ m in equation (|6|) , effectively 
adding up the mass of a completely hydrostatic poly- 
tropic gas, 



M = r\M 



where 



Mdeg 



(7) 



T)M =Vm72 



e deq w d 



eg 



The mass contribution M n d eg due to non-degenerate 
material is given by the difference between equations 
© and ©, 



ndeg 



M - AI, 



4 deg • 



(8) 



We keep track of the time evolution for all three quanti- 
ties in our analysis of white dwarf evolution. Of partic- 
ular interest is the epoch when Mdeg = M n d eg , a time 
that we define as the end of the white dwarf's degener- 
ate phase. 
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2.1.2 Radius 

The radius of the degenerate core is given by 



(9) 



where 



Vr = 



1 



2V2tt \m 



-(3^ 2 ) 2/3 ™- 5/3 ) 1/2 W 



Adding the radial contribution of the non-degenerate 
shell, we find an overall white dwarf radius of 



R(t) 



3M nd 



'■</ 



47T Pq 



R 



dcy 



1/3 



(10) 



Since the contribution to the overall radius due to non- 
degenerate matter is small early on ( only 0.05% for our 
adopted lower bound on density po = 1 g cm" 3 ), an 
approximate form is obtained by substituting £i for ^deg 
in r] r of equation @. This approximation effectively 
assumes that the full, non-truncated density description 
is given by the Lane-Emden equation, so that 



ft = m 



and 



Rdeg 



£,deg 



(ii) 



Because we set a lower bound po on the density at the 
outer surface, the overall radius R is always less than 1Z, 
the radius obtained without assuming a minimum den- 
sity. Calculations of surface temperature in this work 
(see ^2.1.4[) use the radial equation (fTU|) , whereas pre- 
vious calculations (Adams et al. 1998) used the radial 
form given by 

2.1.3 Central Density 

At this point we have expressions for the total mass and 
total radius of the white dwarf, where we have included 
a lower bound po on the density. However, both the 
degenerate mass and radius depend on the central den- 
sity of the white dwarf. The ratio of average density to 
central density of a polytropic system is given by 



£^ 31^)1 

Pc 



6 



(12) 



where here the average density is described in terms of 
the non-truncated mass © and radius HI]), i.e., 



M 



4/3 ttK 3 



(13) 



Note that this average density represents that corre- 
sponding to the entire profile given by the Lane-Emden 
equation, not a profile that is truncated at some non- 
zero density. It is also useful in evaluating equation 
(|13p to invoke the mass-radius relation 

MV? = v 3 nVM c- 3 p: 5 . 

Putting these results together, we find that the cen- 
tral density of the white dwarf can be expressed by the 
following function of time alone, 



c(*) = 



Vp [M 2 (t)G 3 (t)i4(t)] , 



(14) 



where 



Vp 



4^'(a)i; vIvm 



This form provides the correct central density for the 
white dwarf while a degenerate core is present (other- 
wise p c = p a ). 

2.1.4 Luminosity and Temperature 

The luminosity of the white dwarf is provided by proton 
decay and is given by 



L(t) = Tc 2 TM{t) 



(15) 



where J » 2/3 is an efficiency factor {T < 1 due to 
neutrino losses, since such particles do not thcrmalize 
within the star). The surface temperature of the star 
is given by 



T A (t) 



L{t) 
4TTcr sb R 2 



(16) 



where a s b is the Stephan-Boltzmann constant, and R is 
the radius given by equation (TIT))) . Given that R < TZ, 
the effect of including a lower bound on density is to 
provide the star with a somewhat higher temperature, 
effectively shifting its trajectory in the HR diagram to 
the left during late times. 

2.2 Proton Decay 

Proton decay has many possible channels, but the most 
relevant possible processes can be conceptually divided 
into two types. The first type of process takes place 
through unification scale effects, where the intermedi- 
ate particle that drives decay is a GUT scale particle 
that violates baryon number conservation. The second 
type of process is driven by gravitational effects. In 
this case, the simplest description invokes virtual black 
holes, which also violate baryon number conservation, 
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as the effective intermediate particles (for complete- 
ness we note that string theories render this descrip- 
tion overly simplistic). In either case, the estimated 
time scale for proton decay is given approximately by 
the expression 



t P 



(17) 



where r c ~ 10 -31 yr is the light crossing time of the 
proton, mp is the proton mass, and Mx is the mass of 
the intermediate particle that drives the decay process. 

For the case of GUT scale proton decay, Mx ~ 
Mqut ~ 10 16 GeV; in this case, the proton decay time 
scale is about 10 33 yr, roughly the same as the current 
experimental bounds (Super-K 1999). For the case of 
gravitational proton decay, one expects Mx ~ Mpi ~ 
10 19 GeV. This value for the intermediate particle mass 
implies a proton decay time scale of about 10 45 yr, much 
longer than current bounds. Since the suppression of 
baryon number violation at the Planck scale is diffi- 
cult, we expect the proton lifetime to fall in the range 
10 33 yr < tp < 10 45 yr. For the sake of definiteness, we 
take an intermediate time-scale between these limits as 
a working benchmark value so that ~ tp ~ 10 37 yr. 

For the case of gravitationally driven proton de- 
cay with time varying G, an additional complication 
arises. The Planck mass, which sets the effective scale 
for the virtual black holes, scales as Mpi cx G -1 / 2 , so 
the (gravitationally driven) proton decay time scales as 
tp oc G~ 2 . As a result, the proton lifetime is expected 
to increase as the gravitational constant decreases. 

2.3 Chemical Cascade Chain 

This section investigates the variations of fx e = (A/Z), 
which determines, in part, the equation of state. This 
quantity changes with time due to nuclear decay events 
that determine the long term fate of white dwarf stars. 
Initially, a white dwarf star will most likely consist of 
a mixture of heavier elements such as 4 He, 12 C, and 
16 0. On sufficiently long time scales, however, random 
baryon decay events send these particles on a cascading 
path towards 1 H, and eventually total extinction. To 
first order, /Lt e is proportional to the mean value {A/Z) 
of the white dwarf's constituent elements. This ratio 
will be nearly 2 for any white dwarf consisting mostly 
of heavier elements. By the time that about half of the 
original nuclcons have decayed through proton decay 
processes, however, a significant portion of the com- 
posite mass will be in the form of Hydrogen, for which 
fx e = 1. We thus expect fi e to be a decreasing function 
of time, with fj, e — * 1 in the long time limit t — > oo. 



This cascade of particles is categorized below. This 
chain begins with 12 C, and continues down to baryonic 
extinction (Thornton & Rex 2000); note that we could 
start with larger nuclei (e.g. Oxygen), but the latter 
parts of the chain would remain the same. The basic 
reactions include: 
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C-b 



it 



B + 7M 



L B-6^ 10 B + 7M 



9 Bc + lbd 
2 4 Hc+ ^ + 7^ 



3 Be-6 — >2 4 Hc + 7M 



m c -b 



5 He 



3 He 



'-K-b 



X H-^ 



P 
it 



lbd 

2 H + 76d 
2 X H 



lbd 



lbd 



lbd 



(18) 
(19) 

(20) 
(21) 
(22) 

(23) 
(24) 
(25) 



These reaction equations show the nuclei resulting 
from each particular baryon decay event and an en- 
ergy release. Each reaction given above may actually 
include several different reactions en route to the pro- 
duction of the listed stable particles. Among the decay 
products will be some particle with a positive charge, 
often a positron, which annihilates with an electron and 
turns into energy represented by the photons 7^ in the 
reaction equations. Among the baryon decay events 
listed above, (b) signifies either a proton or a neutron 
decaying, (n) signifies neutron decay, and (p + ) signi- 
fies proton decay. The particular daughter particles 
for 10 B and 3 He are dependent upon whether a proton 
or a neutron decays, whereas the other particles listed 
in the chain produce only one variety of products, re- 
gardless of whether a proton or a neutron decays. Al- 
though most reactions listed above have outcomes that 
are independent of which particle decays, the exact path 
taken to achieve the stable particles could vary (and are 
not documented in the chain). 

To model this process, we used a Monte Carlo sim- 
ulation to follow the evolution of the particle species 
population as a function of overall mass for a white 
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Fig. 2 Chemical species mass fraction versus relative 
mass of white dwarf with initial composition of pure 12 C. 
Each curve represents a species of particle within the white 
dwarf during proton decay. The labels indicate the particle 
species, with the exception of 9 Be, which is represented by 
the dashed curve. 

dwarf star undergoing baryonic decay. In these simula- 
tions, both protons and neutrons were allowed to decay 
with equal likelihood; since proton decay (neutron de- 
cay in bound nuclei) has not yet been measured, this 
assumption remains unverified experimentally. If the 
resultant nucleus is unstable to decay, then the decay 
was immediately allowed to take place (since the pro- 
ton decay time-scale is far greater than the nuclear de- 
cay time-scale of any known isotopes of light elements.) 
Figure [3] shows the results of one such Monte Carlo sim- 
ulation for an initial collection of ~ 10 6 atoms of 12 C, 
which were subject to baryon decay (see also AL97). 
The implications of this chemical cascade on the values 
of fx e , and hence the equation of state, are discussed 
below. 

For completeness, we note that the values of [i e can 
be affected by spallation and pycnonuclear reactions. 
In particular, pycnonuclear reactions work to build up 
the abundance of 4 He from smaller elements, effectively 
widening and raising the 4 He curve of Figure [2] at the 
expense of the 3 He, 2 H, and 1 H curves. The effect is 
to increase the value of /j. e and hence the height of the 
curve in Figure [3] near the low mass end. Spallation 
events slowly redistribute protons and neutrons one at 
a time among the nuclei in the star. This redistribution 
happens at a rate that is approximately one per baryon 
decay event. Thus, spallation has a similar effect as py- 
cnonuclear reactions in that it acts to increase the size 
of the nuclei, albeit at a slow rate, and slightly increases 
the value of /j, e near the end of the white dwarf's life. 

Figure 3 shows the ratio function of time, 

where time is given in terms of the white dwarf mass. 



2.2 
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Fig. 3 The ratio fi e as a function of decayed mass fraction 
for white dwarf compositions of initially pure 16 (dashed), 
12 C (solid), and 4 He (dashed). 

These curves are calculated from the Monte Carlo sim- 
ulations described earlier. Initially \x e = 2 for all three 
curves, corresponding to their respective homogenous 
compositions. For a brief period for the 12 C curve, 
this ratio rises above 2 to an approximate value of 2.06 
with nearly 95% of the initial mass remaining. This 
ratio then decreases monotonically towards a value of 1 
(homogeneous 1 H composition). The 16 O curve has a 
much broader hump above the fi e = 2 line and reaches a 
maximum value of approximately 2.07 with nearly 86% 
of the initial mass remaining. These curves would rise 
slightly between the end points by including spallation 
and pycnonuclear effects in the Monte Carlo simulation. 

For simplicity, we use a fitting function that repre- 
sents the \i & curves of Figure [3] with an accuracy of 
~ 1%. The fitting function has the form 

H e (u) = 2 — u + (aiu + d2u 2 + a3M 3 ) exp(<24it) , (26) 

where u = 1 — M/Mq. Although one can find a good 
fit to these curves, there exists no unique set of fitting 
parameters (ai, a.2, a^, 04). Clearly, a\ + a-i + 03 = in 
order to fulfill the requirement that /i e (l) = 1, . But, 
more specifically, for a given £14, there exists no unique 
pair (01, a.2) that provides a fit to within a given accu- 
racy. In Table [T] we provide suggested fitting param- 
eters for three different cases where the initial white 
dwarf composition is taken to be of pure 4 He, 12 C, or 
16 O (note that we give a large number of significant 
figures in this table). 
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Table 1 Fitting parameters for /j, e time evolution curves. 
The Monte Carlo simulations were started with a pure con- 
centration of the chemical species listed in the left column. 
The time evolutions of /x e are recovered by substituting the 
parameters from this table into equation (|26p . 



Species 


CLi 


a-2 


a-3 


0,4 


4 Hc 


-0.952 


3.23 


4.182 


-3.025 


12 C 


3.1218 


-12.9250 


9.24125 


-4.95509 


16Q 


1.8935 


14.3400 


-16.2335 


-5.6545 



3 Results 

Given the model for white dwarf structure constructed 
in the previous section, we consider the long term evo- 
lution of these degenerate stars through the action of 
proton decay and time varying gravitational constant. 
To start, we first consider the case where gravity varies 
but no proton decay takes place. Under these condi- 
tions, the white dwarf grows in radial size with time. 
Both the radius of the degenerate interior and the frac- 
tion of the star that makes up the non-degenerate outer 
layer are increasing functions of time. 

The overall radius increases for two reasons: (1) re- 
ducing gravity's strength effectively reduces the cen- 
tral density of the white dwarf, allowing the degenerate 
core to spread out, and (2) the reduction of the de- 
generate core mass results in an increase in mass of the 
frozen solid lattice (approximately at constant density), 
leading to an increase in its radius. If gravity weakens 
enough, the pressure exerted by gravity at the center 
of the white dwarf decreases below the ~ 10 12 dyne 
cm~ 2 lower bound, completely dissolving the degener- 
ate core into the outer solid lattice. In the long term, 
the white dwarf is left with a radial density that is prac- 
tically uniform and a maximal surface radius. Figure U 
shows the radius as a function of time for this scenario 
with different rates of gravitational time scales t* for 
a IMq, white dwarf with pure 12 C composition. This 
figure shows that under gravitational time variations, 
the white dwarf grows to nearly 100 times its original 
radius before becoming a frozen lattice composed fully 
of non-degenerate matter. 

We also consider the case of constant gravity while 
subjecting the star to baryon decay, which causes the 
star to lose mass. The white dwarf increases in radial 
size until its mass in degenerate matter decreases to 
only about two thirds of the total stellar mass. Af- 
ter this point, the radius and mass of the white dwarf 
shrink until the internal coulomb forces dominate over 
gravity, and the object resembles a rock. This scenario 
is captured by the red curves in Figures and 

Next we consider the effects of both time variations 
in the gravitational constant and baryon decay. The 




t/t. 

Fig. 4 Radius versus time diagram for white dwarf evolu- 
tion with constant 12 C composition (no baryon decay) and 
time varying G with t* = 10 37 yr; p = 0.5 (right), 1.0 (mid- 
dle), and 2.0 (left). The mass is IMq. Radius and time are 
shown in dimensionless form. The dots mark the end of the 
white dwarf's degenerate phase. In all three cases the white 
dwarf grows to almost 100 times its original radius, but the 
time scale on which this growth occurs varies with p. 




Fig. 5 Central density versus time for white dwarf evo- 
lution with initial 12 C composition and = t* = 10 37 
yr; p = 0.0 (red), 0.25 (black), and 0.5 (blue). The start- 
ing mass is IMq. Notice that the central density increases 
before falling off rapidly for all three cases. This feature of 
the central density is due to the initial increase in fi e for 
12 C and thus is expected to occur for an 16 O white dwarf 
as well. Density is shown in g cm -3 and time is given in 
years. 

interplay between these two effects will be the focus 
of the remainder of this section. As shown below, the 
main result is to reduce the ratio of degenerate to non- 
degenerate matter within the star at the time when the 
remnant begins its shrinking phase. 

Figures [6] and show the time behavior of mass and 
radius, respectively, for a 1.0 M Q white dwarf composed 
initially of pure 12 C, subjected to baryonic decay with 
characteristic time scale of T = 10~ 37 yr -1 . The three 
different sets of curves represent various indices of grav- 
itational weakening with time constant t* = 10 37 yr. 
In Figure one can see that the overall mass envelope 
(solid line) for a white dwarf under these conditions is 
independent of the gravitational index p and is fully de- 
scribed by r, the characteristic time for baryon decay. 
The set of red curves show the case of a static gravita- 
tional constant (p = 0). As baryons decay, more mate- 
rial precipitates out of the degenerate core (dashed line) 
to join the outer, non-degenerate layer (dotted line). At 
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time t w 5.4 x 10 37 yr (Ft = 5.4), nearly 99.5% of the 
white dwarf's initial mass has decayed and the mass in 
the degenerate core is less than the mass in the non- 
degenerate outer lattice. This event is marked by the 
solid, red triangles in Figures and [SJ and marks 
the time when the degenerate phase of the star comes to 
an end. This event does not spell the end of the white 
dwarf's life, which goes well beyond t f=a 5.9 x 10 37 yr 
(Tt = 5.9) when the degenerate core disappears entirely 
from the stellar remnant. As shown in Figure this 
time also approximately marks the point when the star 
begins to shrink in radius. At this epoch, R w 1.2 x 10 5 
km, a size that is comparable to that of Uranus or Nep- 
tune. 

Figure [5] shows the H-R diagram, which summarizes 
the evolution of these white dwarf stars. While the star 
remains (mostly) degenerate, the tracks in the H-R di- 
agram are approximately given by L ~ y 12 / 5 . This 
form would be exact for an object that obeys the usual 
white dwarf (degenerate) mass-radius relation and has 
luminosity proportional to its mass. Moderate devia- 
tions from this simple behavior occur due to variations 
in the chemical composition, which affect the equation 
of state (and due to time variations in G - see below). 
After the degenerate phase ends, the tracks in the H-R 
diagram become much steeper with L ~ T 12 . This form 
applies for a constant density object with its luminosity 
proportional to its mass. At this transition point, 99.5% 
of the mass has decayed away, the chemical composi- 
tion for this white dwarf is ~ 99% 1 H, ~ 1% 2 H, with 
traces of 3 Hc and 4 Hc (sec Figure[3]). The remnant now 
amounts to a medium sized planet, composed mostly of 
Hydrogen ice, glowing extremely dim with an output of 
only ~ 2 Watts (Figure [5]) • This huge ball of Hydrogen 
ice fizzles away as the white dwarf continues to fade and 
shrink. Finally, around 2.2 x 10 38 yr (Tt = 22), when 
the mass is ~ 10 21 kg, radius J? = 2.9x 10 2 km, effec- 
tive temperature is T = 1.3 x 10~ 3 K, and luminosity 
L = 1.9 X 10~ 7 Watts, the remnant ceases to be a star. 

This story changes only slightly with the introduc- 
tion of a time varying gravitational constant. The set 
of black curves in these Figures shows what happens 
when the gravitational index is p = 0.25, and blue 
curves represent a gravitational index p = 0.5. For 
the scenario where p = 0.25, the degenerate phase ends 
at t w 4.8 x 10 37 yr (Tt = 4.8) when 99.1% of the 
mass has decayed away. At this time, the star is about 
1.5 x 10 5 km in radius and has a power output of about 
3.3 Watts. This white dwarf is - 98.2% ^ - 1.8% 2 H, 
and ~ 0.07% 3 He by mass, and is thus nearly identi- 
cal chemically to the white dwarf evolving with a static 
gravitational constant. For the scenario where p = 0.5, 
the degenerate phase ends at t f» 4.2 x 10 37 yr (Tt = 4.2) 




Fig. 6 Mass versus time for white dwarf evolution through 
proton decay with an initial 12 C composition and t+ = 
r -i = 10 37 yr . p = o ( red ) 5 o.25 (black), and 0.5 (blue); 
M(t) (solid), Mdeg (dashed), M n d eg (dotted). Mass is shown 
as the remaining fraction and time is shown in years. The 
triangle, square, and circle mark the fractional mass and 
time when the mass in non-degenerate form equals the 
amount of degenerate mass for gravitational index p = 0, 
0.25, and 0.5, respectively. This crossover point defines the 
time when the white dwarf's degenerate phase ends and its 
non-degenerate phase begins. 

when 98.5% of the mass has decayed away. At this 
time, the star is about 1.8 x 10 5 km in radius and has a 
power output of about 5.7 Watts. This white dwarf is 
- 97.1% 1 H, - 2.7% 2 H, and - 0.1% 3 He by mass, and 
is also nearly chemically identical to the white dwarf 
produced using a fixed gravitational constant. 

The white dwarf completes its transformation into a 
rock of hydrogen ice when it becomes optically thin to 
the radiation that originates from within. This condi- 
tion is determined by the star's column density. The 
star contains two types of radiation that must be con- 
sidered when determining this transition. The first type 
is the ~ 235 MeV gamma rays produced in each baryon 
decay event; the stellar remnant becomes optically thin 
to this radiation at an approximate radius of 1 meter 
(for our assumed po). The other type is the radiation 
that comes from the thermalization of these gamma 
rays, which produce photons with wavelength A ~ 10 2 
cm. The white dwarf becomes optically thin to this ra- 
diation at a mass of about 10 21 kg and radius of 2.9 x 10 2 
km. This mass is calculated by using the average den- 
sity of the white dwarf, a value that depends only upon 
coulomb forces between neighboring particles in the Hy- 
drogen lattice and not upon gravity. As a result, the 
transition from white dwarf to rock for the two sce- 
narios with p = 0.25 and p — 0.50 depends only upon 
r, and thus happens at the same time, mass, radius, 
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Fig. 7 Radius versus time for white dwarf evolution 
through proton decay with initial 12 C composition and 
U = r -1 = 10 37 yr; p = (red), 0.25 (black), and 0.5 (blue); 
1Z (dotted), R (solid), Rd eg (dashed). Radius is shown in km 
and time is shown in years. The triangle, square, and circle 
correspond to those shapes in Figure |S] In all cases, these 
degeneracy transitions occur shortly after the white dwarf 
has reached maximum radius and has began to shrink. 




temperature, and luminosity as the case with p = 0. 
Corrections to this approximation could come from ac- 
counting for the exact chemical make up and crystal 
grain structure of the solid lattice, both of which would 
differ slightly for the three indices. 

Figures 151 and [TU1 show curves characterizing the time 
and mass marking the end of the degenerate phase of 
the white dwarf for a continuum of gravitational indices 
spanning 10 orders of magnitude, and seven character- 
istic time scales for gravitational time variations. These 
graphs show the limiting behavior for p -C 1 and p 1 
and for intermediate values. For instance, in the bench- 
mark case of a pure 12 C white dwarf, the dynamical 
range of index parameter is 0.01 < p < 100. If p takes 
on values lower than this range, the fate of the white 
dwarf is dominated by proton decay; if p takes on values 
larger than this range, then the fate of the white dwarf 
is dominated by gravitational time variations. This dy- 
namical range of index parameters broadens and in- 
creases as the product Tt* increases, which indicates 
that gravity must vary rapidly to compete with a pro- 
ton decay. 

Before leaving this section, we briefly consider the is- 
sue of the "starting" mass. First, note that the "start- 
ing" mass in this context is the "final" mass in stel- 
lar evolution calculations. At the end of their nuclear 
burning phase, stars become white dwarfs with masses 
in the range O.O8M < M < 1.4M . To leading order, 
white dwarfs with larger masses must become white 
dwarfs with smaller masses as their constituent baryons 



Fig. 8 H-R Diagram for white dwarf evolution including 
proton decay and time variations in the gravitational con- 
stant. In all cases, the star begins with mass M = 1.0 Mq, 
a pure 12 C composition, and evolves according to t* = 
= 10 37 yr. The three curves show different values of the 
gravitational index p — (red), 0.25 (black), and 0.5 (blue). 
Temperature is shown in Kelvin and luminosity is shown in 
Watts. The triangle, square, and circle correspond to those 
shapes in Figure [6] The dashed line shows the reference 
case of a white dwarf with constant /j, e — 2 and no varia- 
tions in the gravitational constant (p — 0). The lower panel 
shows a close-up of the beginning of the time evolution, for 
the first ~ 10% of the mass loss. During this early phase, 
He > 2, the radius shrinks slightly, and the photospheric 
temperature heats up. 

decay. As a result, white dwarf evolution should not 
depend sensitively on the starting mass (in the long 
term). However, the chemical composition plays a role 
in determining stellar structure, primarily through the 
value of p e , which in turn determines the equation of 
state. White dwarfs with larger initial masses must lose 
a larger percentage of their mass before becoming non- 
degenerate, and hence will experience a greater change 
in chemical composition. One might expect the stellar 
properties at the transition point (from degenerate to 
non-degenerate stars) to depend on the starting stel- 
lar mass. However, the larger white dwarfs begin their 
evolution with larger nuclei (e.g., Carbon and Oxygen), 
whereas small white dwarfs (a few 0.1 Mq) are made 
up mostly of Helium. These two effects largely compen- 
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Fig. 9 Time duration of the degenerate phase of white 
dwarf evolution as a function of the gravitational index p. 
The proton decay time parameter Y = 1CP' !7 yr _1 and the 
gravitational strength time parameter t* = 10 37 yr. The 
initial mass of the white dwarf is 1 Mq . Each set of curves 
contains one black (/x e = 1) and one blue (fi e = 2) dashed 
line and are labeled by an integer i = —3, —2,..., 3 where 
I — log(rti). The time marking the end of the degenerate 
phase is largely independent of p when > and for 
p < 1. For values of p S> 1, the different curves become 
monotonically decreasing, parallel lines. The solid curve 
represents a white dwarf star initially composed of pure 12 C. 

sate, so that white dwarfs of all starting masses become 
non-degenerate with roughly the same properties. 

To illustrate this point, consider white dwarfs sub- 
ject to both baryon decay and time varying G, with the 
standard time scales t* = = 10 37 yr and with p = 
1. We can then examine stellar properties at the time 
when the stars become non-degenerate. For a small 
0.1 M Q white dwarf initially composed of 4 He, this 
transition occurs at t = 1.8 x 10 37 yr (Ft = 1.8)when 
jU e ~ 1.11, the mass M = 1.6 x 10~ 2 M Q , and the ra- 
dius R = 1.9 x 10 6 km. For comparison, for a larger 
1.0 M Q white dwarf initially composed of 12 C, this 
transition occurs at t = 3.3 x 10 37 yr (Ft = 3.3) when 
fj, e « 1.03, the mass M = 3.8 x 10 _2 M Q , and the radius 
R = 2.45 x 10 6 km. These differences are thus at the 
factor of ~ 2 level. 



4 Conclusion 

This paper has constructed a model for the structure 
and evolution of white dwarf stars under the action of 
both proton decay and time variations in the gravita- 
tional constant (section 2). This model includes vari- 
ations in chemical composition (Figure 2 and section 
2.3) which affects the equation of state. The model 
also includes separate accounting for the inner degen- 
erate regions (which obey the polytropic equation of 
state [1]) and the outer non-degenerate regions. 
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Fig. 10 Mass fraction in degenerate matter versus grav- 
itational index p. Each set of curves contains one black 
(fj, e — 1) and one blue (/i e = 2) dashed line and correspond 
to an integer £ = —3, —2,..., 3, from left to right, where 
t = log(rt*). The dynamical range for r = t* has p span- 
ning 4 orders of magnitude centered at p — 1. The solid 
curve represents a white dwarf star initially composed of 
pure 12 C. 

Our main findings can be summarized as follows: 

In the absence of baryon decay, white dwarfs grow 
larger as the strength of gravity decreases (Figure 4). 
The radius reaches a maximum value when the entire 
star becomes non-degenerate; at this time, the radius 
is ~ 100 times larger than its starting value. 

With the inclusion of baryon decay, white dwarfs lose 
mass with time and follow tracks in the H-R diagram 
as shown in Figure [SI In all cases, after an initial tran- 
sient phase, the stars grow dimmer and redder with 
time, and hence move to the lower right in the H-R 
diagram. While the stars remain (primarily) degener- 
ate, the tracks have a shallow slope, with L oc T 12 / 5 . 
After the stars lose enough mass to become primarily 
non-degenerate, the tracks become much steeper, with 
L oc T 12 . These slopes would be exact in the absence 
of changes in the chemical composition; the chemical 
variations and the time variations in G lead to small 
departures from this behavior. 

During the degenerate regime of evolution, the radii 
of these white dwarf stars grow larger due to mass loss 
from proton decay and due to the weakening of grav- 
ity. During the later, non-degenerate phase, the radii 
become smaller. The maximum value of the radius oc- 
curs near the transition between the degenerate and 
non-degenerate regimes, with a value that is typically 
~ 10 times the starting radius (see Figure [7]). 

In addition to understanding the long term fate and 
evolution of white dwarfs, one motivation for this work 
is to understand the larger issue of the future of the uni- 
verse. Previous projections of the future (e.g. AL97, 
Dyson 1979) are predicated on the assumption that 
the laws of physics are known and will not change 
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with time. One important issue is thus the question 
of whether or not the constants of nature change with 
time, and how such variations would change projections 
of our future history. In the case of white dwarf evolu- 
tion, we find that time variation leads to quantitative 
- but not qualitative — changes in the future time- 
line. Time varying G leads to variations in the exact 
mass, size, and composition of white dwarfs as they be- 
come non-degenerate rock-like objects. However, the 
starting states and final states are essentially the same, 
so time variations in gravity do not lead to fundamen- 
tal changes in the overall picture. White dwarfs start 
as degenerate objects with the current value of G, and 
hence are constrained to begin their evolution in nearly 
the same states. At the other end of time, white dwarfs 
cease to act as stars when they become optically thin 
to their internal radiation. At this epoch, the "stars" 
are essentially large rocks of Hydrogen ice, and their 
structure is non-degenerate and largely independent of 
gravity. As a result, time variations in G primarily af- 
fect the intermediate states, and this work shows that 
the modifications are relatively modest (see Figures [5]- 
13). 
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